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ABSTRACT 

We present results from observations of the magnetar IE 1547.0—5408 (SGR J1550— 5418) taken 
with the Chandra X-ray Observatory and the Rossi X-ray Timing Explorer (RXTE) following the 
source’s outbursts in 2008 October and 2009 January. During the time span of the Chandra ob- 
servations, which covers days 4 through 23 and days 2 through 16 after the 2008 and 2009 events, 
respectively, the source spectral shape over the Chandrab&nd remained stable, while the pulsar’s spin- 
down rate in the same span in 2008 increased by a factor of 2.2 as measured by RXTE. This suggests 
decoupling between the source’s spin-down and radiative changes, hence between the spin-down- 
inferred magnetic field strength and that inferred spectrally. The lack of spectral variation during flux 
decay is surprising for models of magnetar outbursts. We also found a strong anti- correlation between 
the phase-averaged flux and the pulsed fraction in the 2008 and 2009 Chandra data, but not in the 
pre-2008 measurements. We discuss these results in the context of the magnetar model. 

Subject headings: pulsars: individual (IE 1547.0—5408, PSR J1550— 5418, SGR J1550— 5418) — stars: 
neutron — X-rays: bursts 

1. INTRODUCTION 

Anomalous X-ray pulsars (AXPs) and soft gamma re- 
peaters (SGRs), though previously thought to be differ- 
ent classes of objects, are now believed to all be strongly 
magnetized neutron stars, known as ‘magnetars’. They 
are characterized by long spin periods of 2-12 s and large 
spin-down rates that imply ultra-strong surface dipole 
magnetic fields of 10 14 — 10 15 G (see reviews by Kaspi 
2007 and Mereghetti 2008) . The most remarkable feature 
of magnetars is their violent outbursts, during which the 
X-ray luminosity could increase by a few orders of magni- 
tude. In the context of the twisted magnetosphere model 
(Thompson et al. 2002), the energy release is due to mag- 
netic field re-arrangement, which may be triggered by 
crustal deformation, and the post-outburst recovery pro- 
vides a powerful diagnostic tool for probing the physical 
conditions of the magnetosphere and the stellar surface. 

Observationally, the transient nature of these events re- 
quires prompt observations; only a handful of studies 
have been carried out with focusing X-ray instruments 
(e.g. Woods et al. 2004; Gavriil et al. 2006; Israel et al. 

2007). In this work, we study recent outbursts from the 
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AXP IE 1547.0— 5408 2 in 2008 and 2009 using observa- 
tions made with the Chandra X-ray Observatory and the 
Rossi X-ray Timing Explorer (RXTE). 

IE 1547.0—5408 was discovered by Lamb & Markert 
(1981) with the Einstein Observatory. Based on its X-ray 
spectrum and infrared flux, Gelfand & Gaensler (2007) 
first suggested the source as a magnetar candidate. The 
detection of radio pulsations by Camilo et al. (2007) di- 
rectly confirmed the pulsar nature of this source; its spin 
period of 2.1s is shorter than any other known mag- 
netars. 3 The pulsar’s spin-down rate as reported by 
Camilo et al. (2007) implies a surface magnetic dipole 
field strength of 2.2 x 10 14 G. Halpern et al. (2008) re- 
ported a high state of the source in 2007 based on XMM- 
Newton observations, and they concluded that an X-ray 
outburst had occurred in between 2006 and 2007. Since 
2008 October 3 (MJD 54742), IE 1547.0-5408 showed 
bursting activity with outbursts detected by Swift (Is- 
rael et al. 2010) and by the Gamma-ray Burst Moni- 
tor (GBM) onboard the Fermi Gamma-ray Space Tele- 
scope (Kaneko et al. 2010). On 2009 January 22 (MJD 
54853), the AXP entered a second active phase. More 
than two hundred bursts were detected within a few 

2 Also known as SGR J1550-5418 or PSR J1550-5418. 

3 See the online magnetar catalog at 
http://www.physics.mcgill.ca/ pulsar /magnetar/ main.html. 
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Fig. 1. — Top: RMS pulsed flux of IE 1547.0—5408 in the 2-10 keV range obtained with RXTE. The red vertical solid lines mark the 
outbursts in 2008 and 2009, and the arrows indicate the Chandra observation epochs. Bottom: total and RMS pulsed Chandra count rates 
in the 1-7 keV band, shown by black squares and red triangles in the upper panels, respectively. The statistical uncertainties are negligible. 
The solid line in 2009 illustrates a power-law fit to the flux decay. Note that the first Chandra observation in 2009 was made with the 
HETG, which precludes a direct comparison of the count rates. Therefore, the data point plotted here is estimated from the spectral 
analysis results. The lower panels show the best-fit blackbody temperature from the PL+BB model, as listed in Table 2. The green dashed 
lines indicate the spin-down rate obtained from the phase- coherent RXTE timing solutions. Uncertainties in v measurements are negligible 
(at 1% level). [See the electronic edition of the paper for a color version of this figure.} 


hours by Swift , INTEGRAL and Fermi GBM (Savchenko 
et al. 2010; Mereghetti et al. 2009; Kaneko et al. 2010). 
Follow-up imaging observations with Swift and XMM- 
Newton revealed dust scattering X-ray rings centered on 
IE 1547.0—5408, from which Tiengo et al. (2010) de- 
duced a source distance of 3.9 kpc. Based on Suzaku 
observations taken 7 days after the 2009 outburst, Enoto 
et al. (2010) reported a hard power- law tail in the spec- 
trum, with a photon index T between 1.33 and 1.55, ex- 
tending up to at least 110 keV. Hard X-ray pulsations 
were also detected by INTEGRAL in the 20-150 keV 
band (Kuiper et al. 2009; den Hartog et al. 2009). The 
pulsed emission has T = 1.55 and the spectral shape re- 
mained stable over the observation period covering from 
day 2 to day 9 after the outburst. 

In this paper, we report on results from Chandra 
and RXTE observations taken after the 2008 and 2009 
events. The observations and results are reported in §2, 
and we discuss their physical implications in §3. We sum- 
marize our findings in §4. 

2. OBSERVATIONS AND RESULTS 
2.1. RXTE Observations and Results 


IE 1547.0—5408 has been monitored regularly with 
RXTE since the 2008 October outburst. Data are col- 
lected with the Proportional Counter Array (PC A) in- 
strument that consists of five collimated xenon/methane 
multi-anode Proportional Counter Units (PCUs). In this 
study, GoodXenonwithPropane mode data were used, 
which give 1-//S time resolution and 256 energy chan- 
nels in the 2-60 keV energy range. We considered events 
from only the top Xenon layer of each PCU to maximize 
the signal-to-noise ratio. Typical exposure times were in 
the range 2-7 ks. In all, we report here on a total of 
55 observations taken between MJDs 54743 and 54911. 
A more complete report on the RXTE data will be pre- 
sented elsewhere (Dib et al, in prep.). 

2.1.1. RXTE Timing 

To determine the pulsar spin parameters, we cleaned 
the RXTE data and selected events between 2-6.5 keV. 
As the source exhibited many short X-ray bursts, we re- 
moved all bursts intervals for this analysis. The pho- 
ton arrival times were first corrected to the solar sys- 
tem barycenter, then binned with 31.25-ms time res- 
olution. The time series were folded at the nominal 
spin period and pulse arrival times were extracted by 
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TABLE 1 

Observation parameters and results 


Obs. 

Date 

MJD 

Obs. ID 

Days since 
outburst 

Exposure 

(ks) 

Count Rate 8, 

(s- 1 ) 

Pulsed Fraction 8, 

2008 October 

1 

2008-10-07 

54746.6 

8811 

4.2 


12.1 

1.35(1) 

0.21(1) 

2 

2008-10-10 

54749.5 

8812 

7.2 


15.1 

1.17(1) 

0.22(1) 

3 

2008-10-18 

54757.6 

8813 

15.2 


10.1 

1.15(1) 

0.33(1) 

4 

2008-10-21 

54760.8 

10792 

18.4 


10.1 

1.07(1) 

0.35(1) 

5 

2008-10-26 

54765.1 

8814 

22.8 


23.0 

0.99(1) 

0.31(1) 

2009 January 

6 

2009-01-23 

54855.0 

10185 

2.0 


10.1 

0.95(l) b 

0.09(1) 

7 

2009-01-25 

54856.7 

10186 

3.7 


12.1 

3.1(2) 

0.09(1) 

8 

2009-01-29 

54860.8 

10187 

7.8 


13.1 

2.5(1) 

0.13(1) 

9 

2009-02-06 

54868.6 

10188 

15.6 


14.3 

2.2(1) 

0.12(1) 


a In the 1 — 7keV energy range. 

k Obs. 6 was made with the HETG, which results in a reduced count rate. 


cross-correlating with a template. These arrival times 
were then fitted to a simple phase- coherent timing model 
using the TEMPO software package 4 . Details on the 
phase folding and periodicity search techniques are de- 
scribed in Dib et al. (2009). Our best-fit timing solution 
gives a spin frequency v = 0.48277818(5) Hz, frequency 
derivative v = —6.6(1) x 10 -12 s -2 and second frequency 
derivative v = —6.2(1) x 10 _18 s“ 3 covering the 2008 
Chandra epochs (reference epoch of MJD 54743.0), and 
v — 0.4825962(4) Hz, v = —5.17(5) x 10 -12 s -2 with neg- 
ligible v for the 2009 Chandra epochs (reference epoch of 
MJD 54854.0). These values are consistent with those 
reported by Israel et al. (2010) and Kaneko et al. (2010) 
based on the 2008 Swift and 2009 Fermi observations, 
respectively. 

2 . 1 . 2 . RXTE Pulsed Fluxes 

Once the timing solution was obtained, we re- 
generated pulse profiles in the 2-10 keV band, and cal- 
culated the RMS pulsed flux according to the formula in 
Dib et al. (2008), using 7 harmonics. The results are plot- 
ted in Figure 1 and reveal a complicated flux evolution. 
In this study, we focus mainly on time periods near the 
Chandra epochs, and a full analysis of the RXTE data 
will be presented by Dib et al. (in prep.). Approximately 
11 days after the initial 2008 October 3 (MJD 54742) 
trigger, during which the pulsed flux decayed roughly 
monotonically, the pulsed emission abruptly increased 
again by 80% between two RXTE observations taken 
on MJDs 54751.2 and 54752.1, then decayed monotoni- 
cally until it reached a minimum value on MJD 54786.2. 
Since the onset of this second event is not resolved, we 
report an upper limit of lday for the rise time. This 
second flux enhancement decayed monotonically for ~34 
days; the pulsed count rate in this period can be param- 
eterized by an exponential fall-off with 1/e decay time of 
~ 25 days although a linear decay is also consistent with 
the data. Interestingly, the 2009 event is far less dra- 
matic as seen by RXTE , with a relatively small increase 
in the pulsed flux observed. 

2.2. Chandra Observations and Results 

The 2008 outburst of IE 1547.0— 5408 triggered a series 
of Chandra observations through our Target of Oppor- 

4 http://www.atnf.csiro.au/research/pulsar/tempo. 


tunity (ToO) program. Five total pointings were made 
on days 4, 7, 15, 18 and 23 after the outburst, with ex- 
posures ranging from 12 to 23 ks. Data were taken with 
the ACIS-S detector in continuous clocking (CC) mode, 
which has a time resolution of 2.85 ms. For the 2009 
event, another Chandra ToO program followed the out- 
burst, with four ACIS-S CC mode observations taken on 
days 2, 4, 8 and 16. The first exposure in 2009 was taken 
with the High Energy Transmission Grating (HETG) . In 
this study, we include only the zeroth-order events in this 
latter data set. A summary of observation parameters is 
provided in Table 1. 

We carried out the Chandra data reduction using 
CIAO 4.2 with CALDB 4.2.0. Source counts in all the 
data sets were extracted from a 3"-diameter aperture 
and background counts were from the whole chip ex- 
cluding the central 40" region (i.e. a total width of 7'.7). 
We note that although there is a nearby source XMMU 
J155053. 7— 541925 at 65" away from the pulsar (Gelfand 
& Gaensler 2007), the projected separation is greater 
than 50" in all our CC-mode data. Therefore, it does 
not contaminate the source counts. The source count 
rates in the 1-7 keV energy range are reported in Table 1 
and plotted in Figure 1, in which the backgrounds have 
been accounted for, although they are less than 0.5%. 
The count rates are well below the telemetry limit of 
Chandra , and pileup is negligible due to the short frame 
time of the CC-mode exposures. While the source flux 
changes between epochs, we found no short-term vari- 
ability within any individual exposure. Employing a test 
algorithm suggested by Gregory & Loredo (1992), we 
found < 10% probability of source variability within in- 
dividual observations, which is not statistically signifi- 
cant. The flux decay in 2009 can be modeled by a power 
law of index a = — 0.21d=0.01, which is plotted in Fig- 
ure 1. However, such a simple relation is not observed 
in the 2008 flux evolution. In particular, the count rates 
are very similar between the second and third exposures, 
which is likely related to the enhanced pulsed flux de- 
tected by RXTE. 

2 . 2 . 1 . Chandra Timing 

We applied a barycenter correction to the Chandra 
data and then folded the photon arrival times accord- 
ing to the ephemeris above. The resulting light curves 
are shown in Figure 2. The pulse profiles in 2008 show 
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2008 2009 



Fig. 2. — Pulse profiles of IE 1547.0—5408 in 1-7 keV obtained 
from the Chandra observations, using 64 phase bins. The RMS 
pulsed fractions from Table 1 are indicated. 

some hints of a multi-peak morphology in the first two 
observations, then evolve into a single peak. By contrast, 
the 2009 profiles exhibit a broad peak at first, with a sec- 
ond peak emerging by the third observation. We found 
no obvious energy dependence of the pulse shape across 
the Chandra 0.5-10 keV band. A direct comparison be- 
tween the profiles in 2008 and 2009 indicates a much 
higher pulse modulation in the former. We estimated 
the RMS pulsed fractions (PFs) in 1-7 keV, and the re- 
sults are listed in Table 1. We observe a clear trend of 
increasing pulse modulation as the source recovers after 
the outbursts. 

To look for energy dependence of the modulation, we 
estimated the PFs in the soft (1-3 keV) and hard (3- 
7keV) energy bands separately. In the 2008 data, the 
latter show a systematically higher PF, with a difference 
ranging from APF=0.04-0.09 (i.e. a 20-30% change), 
which is statistically significant given the ~0.01 uncer- 
tainties. However, such an energy dependence is not ob- 
served in 2009, with the PFs in the two bands being 
consistent with each other. 

2.2.2. Chandra Spectroscopy 

The Chandra spectra of IE 1547.0—5408 were ex- 
tracted using the tool psextract in CIAO, then binned 
such that every bin has a S/N of at least 10. We 
performed the spectral fits in the 0.5-7 keV range with 
XSPEC vl2.6.0. All nine datasets were fitted jointly with 
a single absorption column density (AT H ). We also tried 
fitting different 7V H values for the 2008 and 2009 data, 
and confirmed that they are consistent. We started with 
simple models including an absorbed blackbody (BB) 
and an absorbed power-law (PL), but obtained very poor 
fits (reduced x 2 values over 1.5). An absorbed black- 


2008 BB+PL model 



2009 BB+PL model 



Fig. 3. — Best-fit blackbody plus power-law model fit to the 2008 
and 2009 Chandra spectra (upper panels) with the corresponding 
residuals (lower panels). Different observations are shown by dif- 
ferent colors. The corresponding best-fit spectral parameters are 
listed in Table 2. [^ee the electronic edition of the paper for a color 
version of this figure.] 

body plus power-law (BB+PL) model gives much better 
fits and the results are listed in Table 2 and plotted in 
Figure 3. As shown in the Figure, the fit residuals sug- 
gest a hint of a spectral feature ~ 1.3 keV, which is more 
obvious in 2009 than in 2008. However, the significance 
is only ~1 a and deeper exposures are needed to confirm 
this. 

In addition to the BB+PL model, we also considered 
more physical models that account for the Compton up- 
scattering of thermal photons in the magnetosphere. We 
tried fitting the Resonant Cyclotron Scattering (RCS; 
Lyutikov & Gavriil 2006; Rea et al. 2008) and the Sur- 
face Thermal Emission and Magnetospheric Scattering 
(STEMS; Ozel 2003; Giiver et al. 2006) models to the 
data. In the latter, the gravitational redshift is fixed at 
z = 0.306 during the fit, corresponding to the canon- 
ical neutron star mass of 1.4 M 0 and radius of 10 km. 
While these models fit the 2008 data reasonably well, 
the 2009 spectra clearly require an additional hard com- 
ponent. Therefore, we added a PL to the spectral model 
in 2009, with T fixed at 1.33 according to the Suzaku 
results (Enoto et al. 2010). 5 Comparing to the BB+PL 
fit above, these models provide a similar goodness- of- fit 

5 The INTEGRAL results also suggest that the hard-band PL 
spectral index remained stable over the period of the Chandra ob- 
servations (den Hartog et al. 2009). 
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TABLE 2 

Phase- averaged spectral parameters of IE 1547.0—5408 for different models 


Obs. 

Ah 

B 

kT 

R 

r 

jabs a 

F p i/F th b 

Xt/v 


(10 22 cm -2 ) 

(10 14 G) 

(keV) 

(km) 


(10 -11 ergs 
S -1 cm -2 ) 



Blackbody + Power-law model 

2008 









1 

4.1 ±0.1° 


0.58 zb 0.03 

2.0 ±0.3 

9 0+0.3 
Z- ^— 0.4 

1.92 ±0.03 

1 4 +1 -° 
1 - 0.6 

1.08/1342° 

2 

4.1 ±0.1° 


0.57 zb 0.02 

9 -I +0.3 
Zi - 0.2 

9 0+0.3 
z -' 3 — 0.4 

1.64 ±0.02 

x x+ 0 - 8 

1.08/1342° 

3 

4.1 ±0.1° 


0.60 zb 0.04 

X 7 +° - 4 

A- —0.3 

2.8 ±0.3 

1.57 ±0.03 

2 . 6 +X? 

1.08/1342° 

4 

4.1 ± 0.1 c 


0 e;e;+ 0- 03 
U.DD_o .02 

9 9 + 0.4 
ZZ -0.3 

2 4 +°- 3 
Z - -0.5 

1.55 ±0.03 

1 2+ 1 - 1 
1 -Z — 0.6 

1.08/1342° 

5 

4.1 ± 0.1 c 


0.57 zb 0.02 

1.8 ±0.2 

9 O+ 0.2 

1.34 ±0.02 

2 X - *" 1 ' 2 
Z - 0.8 

1.08/1342° 

2009 







6 

4.1 ± 0.1 c 


0.49 zb 0.03 

4 4+0-8 
-0.7 

2 Q+0.3 
zu -0.4 

5.11 ±0.09 

1 6 +1 ' 8 
AD - 0.8 

1.08/1342° 

7 

4.1 ± 0.1 c 


0.46 zb 0.02 

4 7+0-7 
^■'-0-5 

2.1 ± 0.2 

4.61 ± 0.05 

1 8 +1 ‘° 
iO - 0.6 

1.08/1342° 

8 

4.1 ± 0.1 c 


0.47 zb 0.02 

4 0+ 0 - 6 
^■ U -0.5 

2.3 ±0.2 

3.67 ±0.04 

9 9+1-3 
Z-Z — 0.8 

1.08/1342° 

9 

4.1 ± 0.1 c 


0 44 + 0 02 
-0.03 

4 1 + 0 - 8 
^■ A - 0.6 

2.4 ±0.2 

3.22 ± 0.04 

9 7 +I -8 
z - ‘ - 1.0 

1.08/1342° 

RCS + hard Power-law model 

2008 









1 

4.0 ± 0.1 c 


0 

1 + 
o o 
o b 

4^ Or 



1.91 ±0.03 


1.09/1338° 

2 

4.0 db 0.1 c 


0.47 ±0.04 



1.63 ±0.02 


1.09/1338° 

3 

4.0 ± 0.1 c 


U.OO — o 05 



1.56 ±0.03 


1.09/1338° 

4 

4.0 db 0.1 c 


0 43 + 0 07 
u-^o_ 0 .03 



1.53 ±0.03 


1.09/1338° 

5 

4.0±0.1 c 


0.41 ±0.04 



1.33 ±0.02 


1.09/1338° 

2009 








6 

4.0 ± 0.1° 


0.32±°;°s 


1.33 d 

5.11 ±0.09 

0.37 ±0.08 

1.09/1338° 

7 

4.0 db 0.1 c 


O-i+o.i® 


1.33 d 

4.62 ± 0.05 

°- 3 + 0 .i 2 

1.09/1338° 

8 

4.0 =L 0.1 c 


°-32t° 0 + 


1.33 d 

3.67 ±0.03 

°- 30 io.09 

1.09/1338° 

9 

4.0 ± 0.1° 


0-13t° 0 + 


1.33 d 

3.23 ±0.04 

°- 27 io.08 

1.09/1338° 

STEMS + hard Power-law model 

2008 









1 

4.5 =L 0.1 c 

2 fil + 0 - 10 
z - o± -0.08 

0 Q1 Q+ 0 011 

U.0±0_Q QQg 

10 d 


1.90 ±0.03 


1.08/1329° 

2 

4.5 d= 0.1° 

9 HQ+O-IO 
z - oy -0.07 

0 314+ 0011 
U. 0.005 

10 d 


1.62 ± 0.02 


1.08/1330° 

3 

4.5 d= 0.1° 

2.57 zb 0.06 

n 9 C+ 0-02 

U.Z»_0.06 

10 d 


1.56 ±0.03 


1.08/1331° 

4 

4.5 ± 0.1° 

o 4+0 - 2 

*- 4 -o.3 

0 Qc;+ 0 - 05 

U.OO_o 02 

10 d 


1.53 ±0.03 


1.08/1332° 

5 

4.5 zb 0.1° 

9 ft 9 + 0-06 

z.uz— o 04 

n qn+0.005 

U.OU _ 0 07 

10 d 


1.33 ±0.02 


1.08/1333° 

2009 









6 

4.5 zb 0.1 c 

5-4 ±?:I 

0 46 +0 09 

10 d 

1.33 d 

5.08 ±0.09 

0 21 +0 09 
UZ 1 - 0.11 

1.08/1334° 

7 

4.5 zb 0.1° 

3-01 + o'o9 

0.306 ± 0.003 

10 d 

1.33 d 

4.60 ±0.05 

0.16 ± 0.02 

1.08/1335° 

8 

4.5 zb 0.1 c 

3 fi+ 0 - 2 

n qi 9 + 0-012 
U.olz_Q 002 

10 d 

1.33 d 

3.66 ±0.04 

0 X 2 +0 ' 03 
UiZ - 0.02 

1.08/1336° 

9 

4.5 zb 0.1° 

o q+0.6 

0 so - * -0 02 

U .OU_ 0 oil 

10 d 

1.33 d 

3.23 ±0.04 

0 1 c;+ 0 - 02 
Uitl - 0.01 

1.08/1337° 


Note. — All uncertainties quoted axe 90% confidence intervals (i.e. 1.6<r). 
a Absorbed flux in 0.5-7 keV range. 

k Unabsorbed flux ratio between the power-law and the primary (i.e. blackbody, STEMS, or RCS) components in 0.5-7 keV 
range. 

c All nine datasets were fitted jointly with a single column density. 

^ Held fixed in the fit. 


in terms of the reduced x 2 values. Table 2 lists the key 
parameters of the best-fit models. The scattering optical 
depth r is around 1-2 in 2008 and > 3 in 2009, and the 
thermal velocity /3 of the electrons is ~0.4-0.5 in 2008 
and ~ 0.2 in 2009 (see Lyutikov & Gavriil 2006, for a 
detailed definition of these parameters). 

3. DISCUSSION 

In this paper, we have reported on Chandra observa- 
tions of IE 1547.0—5408 immediately following its 2008 
and 2009 outbursts, along with RXTE timing and pulsed 
flux behavior following the 2008 outburst and through- 
out the 2009 event. Next we discuss these observations 
in the context of the magnetar model. 

3.1. Spectral and Spin Evolution 


In the twisted magnetosphere model of magnet ars 
(Thompson et al. 2002), the observed X-ray luminosity 
of a magnetar is determined both by its surface temper- 
ature and by magnet ospheric currents, the latter due to 
the twisted dipolar field structure. The surface temper- 
ature in turn is determined by the energy output from 
within the star due to magnetic field decay, as well as 
on the nature of the atmosphere and the stellar mag- 
netic field strength. This surface thermal emission is 
resonantly scattered by the current particles, thus re- 
sulting in an overall spectrum similar to a Comptonized 
blackbody (e.g. Lyutikov & Gavriil 2006; Rea et al. 2008; 
Zane et al. 2009). In this model, the greater the twist an- 
gle, the greater the scattering, the harder the spectrum, 
and the greater the X-ray luminosity Lx ■ In addition, 
the surface heating by return currents is believed to con- 
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tribute substantially to Lx, at least at the same level as 
the thermal component induced from the interior field 
decay (Thompson et al. 2002). Magnetar outbursts in 
this picture occur with sudden increases in twist angle, 
consistent with the generic hardening of magnetar spec- 
tra during outbursts (e.g. Kaspi et al. 2003; Woods et al. 
2004; Israel et al. 2007). 

Other observational evidence provided in support 
of the twisted magnetosphere model as proposed by 
Thompson et al. (2002) is a correlation between mag- 
netar spectral hardness and spin- down- inferred mag- 
netic field strength B, when comparing different sources 
(Marsden & White 2001; Kaspi & Boydstun 2010). In 
this case, B is an observational proxy for the magne- 
tospheric twist angle. On the other hand, some mag- 
net ars have shown dramatic spin-down rate variations, 
with order-of-magnitude changes in v seen on a variety 
of time scales (e.g. Gavriil & Kaspi 2004; Woods et al. 
2004). The origin of these variations is unknown. Nev- 
ertheless, in the context of the twisted magnetosphere 
model, a varying twist angle might naively be expected 
to be accompanied by a changing v (due to changing ef- 
fective B ) , and corresponding spectral and flux changes. 
On the other hand, some decoupling between i > and the 
radiative behavior might be expected, particularly as the 
spin-down is affected most by a narrow field-line bun- 
dle near the light cylinder, whereas the radiation orig- 
inates from the surface. Field- line twists likely propa- 
gate outward (Thompson et al. 2002), which suggests 
the radiative changes should occur prior to v changes 
(Beloborodov & Thompson 2007). 

In contrast to the picture in which a magnetar outburst 
is accompanied by an enhanced magnetospheric twist, 
Ozel & Giiver (2007) suggest that in outburst, the mag- 
netosphere may be stable, with radiative evolution being 
due to changes in the surface thermal emission. Using 
a spectral model consisting of a resonant Comptonized 
atmosphere- modified blackbody (the STEMS model), 
fits to data for XTE J1810— 197 (Giiver et al. 2007) result 
in the spectrally inferred B being stable, with all radia- 
tive changes being due to changes in the surface thermal 
emission. 

For IE 1547.0—5408, the source spectrum only showed 
significant variability over a short period of time ( rsj 1 day) 
after the 2008 and 2009 outbursts (Israel et al. 2010; 
Scholz et al. in prep.), but remained stable over the 
Chandra observations, during which the flux changed 
greatly (see Table 2). The latter seems opposite to the 
hardness/flux correlation predicted by the twisted mag- 
netosphere model. More intriguing are the timing re- 
sults. As described in §2.1.1, we found that the spin- 
down rate |z>| increased by a factor of 2.2 between the 
first and last Chandra observations in 2008 (see §2.1.1), 
a substantial change even by magnetar standards. For 
purely dipole spin-down, this naively implies a ~50% in- 
crease in the effective surface dipole field strength B, 
from an initial spin- inferred value of 2.8 x 10 14 G at the 
epoch of the first Chandra observation, to a value of 
4.1 x 10 14 G at the epoch of the last. This is contrary to 
SGR 1806—20, in which the spectral response lagged be- 
hind the torque variation, suggesting some hysteresis in 
the system (Woods et al. 2007). In our case, the lack of 
associated spectral changes is unexpected in the twisted 


magnetosphere model, unless the spin-down is decoupled 
from the site of the radiative events, as suggested by 
Gavriil & Kaspi (2004) for IE 1048.1—5937. Our results 
call into question the reliability of spin- down- inferred B 
values being used when comparing those inferred spec- 
trally. 

Indeed, for the 2008 Chandra observations of 
IE 1547.0—5408, spectral fits using the STEMS model 
yield a very similar B-field value for all the observations 
(see Table 2) in spite of the strongly varying z>. One 
way to reconcile this is to interpret the 5-field measured 
spectroscopically as higher-order multipoles in localized 
X-ray-emitting regions rather than the global dipole field 
responsible for spin-down. However, this is hard to rec- 
oncile with the fact that the STEMS B values are lower 
than the spin-inferred value at the epoch of the last 

2008 Chandra observation. As an alternative, the extra 
spin-down torque could be attributed to particle winds 
(Harding et al. 1999). Comparing the spin-down rates of 
IE 1547.0-5408 in 2007 (Camilo et al. 2008), 2008 and 

2009 (§2.1.1), it is obvious that the spin-down torque 
changed drastically between these epochs. Torque varia- 
tions have been observed in magnetars and ordinary ra- 
dio pulsars (e.g. Kaspi et al. 2003; Gavriil & Kaspi 2004; 
Kramer et al. 2006; Lyne et al. 2010; Livingstone et al. 
2010). These could be related to changes in plasma con- 
ditions in the magnetosphere, which may not necessarily 
have any observable effects in the radiative properties 
(see Livingstone et al. 2010). Based on Equation 9 in 
Harding et al. (1999), a factor of 2.2 increase in |z>|, and 
hence in 5, requires a steady wind luminosity ( L p ) of 
1.5 times larger than the dipole spin-down luminosity, 
implying L p « 1.5 x 10 35 ergss -1 . For a typical X-ray 
efficiency of below 1%, the particle- induced flux would be 
below 10~ 12 ergs cm -2 s -1 , hence is likely undetectable, 
given the source’s much higher flux (see Table 2). 

In any case, the absence of spectral variations in the 
presence of flux changes remains puzzling whether or not 
the magnetospheric twist angle varied in the outburst. 
Moreover, a reasonable model will need to explain why z> 
changed drastically following the 2008 event, but stayed 
constant in 2009. This may reinforce the requirement of 
decoupling between the spin-down and the source of the 
radiative changes, hence presumably between the spin- 
down-inferred B and that inferred spectrally. 

3.2. Flux Evolution During the 2008 and 2009 Events 

The 2008 and 2009 events exhibited very different flux 
evolutions. Immediately after the 2009 outburst, the 
persistent flux increased by a factor of ~ 500 (Scholz 
et al. in prep.), while the pulsed flux evidently showed 
only very little variation (less than a factor of 2). We 
found a monotonic flux decay during the recovery, with 
a power-law index of —0.21 ± 0.01, which is compa- 
rable to -0.306 ± 0.005 for CXOU J164710.2-455216 
(Woods et al. 2010), but not as steep as —0.69 ± 0.03 for 
IE 2259+586 or -0.92 ± 0.02 for XTE J1810-197 (Zhu 
et al. 2008; Woods et al. 2005). 

In contrast, the 2008 event is less energetic; the total 
flux increased by a factor of ~100 (Israel et al. 2010), 
but our RXTE results reveal a pulsed flux variation of 
a factor of ~4, far greater than in the 2009 event. Also, 
the flux decay in 2008 showed a more complicated evolu- 
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tion. As is clear in Figure 1, we observed an additional 
flux enhancement around MJD 54752, ~11 days after the 
initial trigger, 6 lasting for ~30 days. 

Eichler & Shaisultanov (2010) suggest that radiative 
outbursts in magnetars could generally be preceded by 
glitches, with the delay between the two events depend- 
ing on the depth at which the glitch- induced energy re- 
lease occurs. For the initial event in 2008, we are un- 
able to tell whether a glitch preceded the radiative out- 
burst, due to the lack of RXTE observations prior to the 
outburst. However, we can rule out any glitch between 
the initial 2008 event and the second flux enhancement 
10 days later. It is possible that the initial event actu- 
ally involved glitches occurring in two different places in 
the stellar interior, at substantially different depths, such 
that the delays from the glitch to the X-ray enhancement 
were different, but we note that this picture does not ex- 
plain the sharp rise of the second flux enhancement. 



Abs. flux (ergs s _1 cm -8 ) 

Fig. 4. — RMS pulsed fraction vs. absorbed flux of 
IE 1547.0—5408 in the 1-7 keV band. The filled triangles and 
squares show the 2008 and 2009 Chandra observations, respec- 
tively. The green open circle and red dots represent results 
from XMM-Newton and Swift , respectively (Halpern et al. 2008; 
Olausen et al. in prep.; Scholz et al. in prep.), with the XMM- 
Newton epochs labeled. Note that the 2007 and 2008 XMM- 
Newton data points are only upper limits, since the values reported 
in the literature are area pulsed fraction estimates, which are larger 
than the RMS pulsed fractions by a constant which depends on the 
pulse shape. The dashed line shows a power-law fit to the Chandra 
data. [See the electronic edition of the paper for a color version of 
this figure.] 


3.3. Pulsed Fraction Evolution 

Our results in §2 clearly indicate a strong anti- 
correlation between the pulsed fraction and the phase- 
averaged X-ray flux, at least during the 2008 and 2009 
outbursts. This is plotted in Figure 4, and suggests an 
approximate power-law relation between the two observ- 
ables. The trend is also supported by the XMM-Newton 
and Swift measurements taken in the same period 
(Olausen et al. in prep.; Scholz et al. in prep.). Similar 

6 Although this second enhancement is not reported by Israel 
et al. (2010) in their study of the Swift observations covering the 
same period, we have re-analysed the same Swift dataset and con- 
firmed the pulsed flux enhancement we observe with RXTE and 
Chandra (Scholz et al., in prep.). 


anti-correlations have been observed in IE 1048.1—5937 
and CXOU J164710.2-455216 (Tam et al. 2008; Israel 
et al. 2007), while positive correlations were found in 
XTE J1810— 197 and IE 2259+586 (Gotthelf et al. 2004; 
Zhu et al. 2008). This variety of behaviors is consistent 
with the picture in which, from source to source, the lo- 
cation and geometry of the region on the star affected 
in the outburst are different. Previous studies proposed 
that an anti-correlation could be the consequence of an 
increased emitting area due to an outburst, such that 
part of the hot spot becomes visible at any phase, thus 
reducing the pulse modulation (e.g. Halpern et al. 2008). 
We note that this scenario depends critically on the loca- 
tion of emission zone on the stellar surface as well as the 
viewing geometry (see Bogdanov et al. 2008); it may be 
possible to obtain either a monotonic increase or decrease 
of pulsed fraction for the same area ‘hot spot’ depending 
on its location on the stellar surface. 

We point out that the trend of decreasing pulsed frac- 
tion for increasing phase- averaged X-ray flux observed 
in the Chandra data does not seem to hold for other 
observations of IE 1547.0—5408 before 2008. Based on 
XMM-Newton exposures, Halpern et al. (2008) reported 
area PFs of 15% in quiescence in 2006 and 7% in the high 
state in 2007. As shown in Figure 4, these values deviate 
significantly from the 2008-2009 trend. The discrepancy 
seems too large to be reconciled by a difference in the 
instrument response, and is even larger if the differing 
estimate methods are accounted for. It is possible that 
the outburst in 2008 induced some permanent changes in 
the field configuration or emission geometry. The pulse 
profiles shown in Halpern et al. (2008) also appear to 
have a different shape from the ones shown in Figure 2, 
providing further support to this picture. 

4. CONCLUSIONS 

We have presented results from Chandra and RXTE 
observations of IE 1547.0—5408 following its 2008 and 
2009 outbursts. These allow a direct comparison between 
the two events. We found that over the 2008 Chandra 
observation epochs, the pulsar spin-down rate increased 
by a factor 2.2, in the absence of corresponding spec- 
tral changes, whereas such variation in v is not observed 
after the more energetic 2009 event. This provides evi- 
dence of decoupling between magnetar spin and radiative 
properties. The absence of spectral changes simultane- 
ous with significant flux decay is surprising for models of 
magnetar outbursts. Our results also revealed a strong 
anti-correlation between the pulsed fraction and phase- 
averaged flux of the source. While both 2008 and 2009 
data follow the same trend, pre-2008 measurements show 
significant deviation, suggesting that the 2008 outburst 
may have induced permanent changes in the emission 
geometry. Finally, we note that IE 1547.0—5408 demon- 
strated significant spectral changes only within the first 
day after the 2008 and 2009 events, which highlights the 
importance of prompt observations in future studies for 
understanding post-outburst relaxation of magnetars. 
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